Prediction  of  Atmospheric  MTF  and 
Application  to  Image  Restoration,  Based 
on  Meteorological  Data 


3rd  Quarterly  Report 

(May  1994) 


N.  S.  Kopeika,  Principal  Investigator 
Department  of  Electrical  and  Computer  Engineering 
Ben-Gurion  University  of  the  Negev 
Beer-Sheva,  Israel 


REPORT  DOCUMENTATION  PAGE 


Form  Approved  0MB  No.  0704-01 88 


Public  reporting  burden  for  «i.  coHedkx.  of  intormati^i.  esw^ed  to 

gMharing  «.d  mainbrining  the  date  needed,  and  eomptet|ng  and  r^^g  me  Opetatlona  and  Reports,  1215J^erson 

collection  of  Information  inciudina  fugoestwns  for  raducinfl  this  burden  to  WashmQton  Headquarters  Services,  D  /r\7Aii_niftR\  wa^hinntnn  np  20503 

22202^302.  and  to  the  Office  of  Management  and  Budget.  PaperworK  Reduction  Proiect  (0704-0188).  Washington.  DC  20503. - 

....  I  H  r^ATr:  *  To  DCD/^DT  TVPP  AND  DATES  COVERED 


1 .  AGENCY  USE  ONLY  (Leave  blank) 


2.  REPORT  DATE 

Nay  1994 


3.  REPORT  TYPE  AND  DATES  COVERED 
Final  Report 


4.  TITLE  AND  SUBTITLE 

Prediction  of  Atmospheric  MTF  and  Application  to  Image 
Restoration,  Based  on  Meteorological  Data. 

3rd  Quarterly  Report  _ _ _ _ _ _ — 


6.  AUTHOR(S) 


I.  Dror  and  N.S.  Kopeika 


5.  FUNDING  NUMBERS 

F6170893W0896 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 


Ben  Gurion  University  of  the  Negev 
Department  of  Electrical  and  Computer  Engineering 
Beer  Sheva,  Israel  84105 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

SPC-93-40603 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 


EOARD 

PSC  802  BOX  14 
FPO  09499-0200 


10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 

SPC-93-4060 


11.  SUPPLEMENTARY  NOTES 


12a.  DISTRIBUTION/AVAILABILITY  STATEMENT 


Approved  for  public  release:  distribution  is  unlimited. 


12b.  DISTRIBUTION  CODE 
A 


13.  ABSTRACT  (Maximum  200  words) 

Predictions  of  atmospheric  transmittance  in  desert  aerosol  environments  using  MODTRAN  code  diverge  significantly  from 
prediction  of  the  desert  particulate  size  distribution  is  required  in  order  to  predict  atmospheric  scattenng  and  absorption  parameters  It  is  also  ess^iai 
to  the  prediction  of  the  aerosol  atmospheric  modulation  transfer  function  which  is  often  the  dominant  component  of  the  overall  atmosphenc  l^F^ 
Recently  an  effort  to  predict  statistics  but  not  size  distribution  according  to  simple  weather  parameters  has  been  made  for  coarse  desert  aerosols, 
quantitative  analysis  of  the  desert  particulate  size  distribution  models  was  also  performed.  In  this  research.  ^^® 

by  optical  counters  are  related  to  weather  parameters.  Known  statistical  and  analytical  models  such  as  MODTI^N  relate  the  size 
parameters  only  to  relative  humidity  for  continental  atmospheres.  Although  humidity  has  a  significant  role  in  the  prediction  of  aerosol  size  statistics 
other  weather  parameters  are  seen  here  to  strongly  influence  also  the  size  distribution  parameters.  Compansons  such  as  the  above  can  be  used  to 
predict  under  which  conditions  the  MODTRAN  aerosol  models  have  good  or  poor  accuracy.  It  is  also  hoped  that  they  will  ead  to  '"[“f 
MODTRAN,  improving  the  accuracy  of  the  humidity  dependence  as  well  as  by  incorporating  other  meteorological  parameters  into  the  MOD 

prediction  models.  ^ . 


14.  SUBJECT  TERMS 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 

UNCLASSIFIED 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

UNCLASSIFIED  _ 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 

UNCLASSIFIED _ 


15.  NUMBER  OF  PAGES 


16.  PRICE  CODE 


20.  LIMITATION  OF  ABSTRACT 


UL 


NSN  7540-01-280-5500 


Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std.  239-18 


A  n'i 


Prediction  of  Atmospheric  MTF  and 
Application  to  Image  Restoration,  Based 
on  Meteorological  Data 

Third  Quarterly  Report 

This  quarter  work  continued  once  again  on  image  restoration,  and  on  aerosol 
size  distribution  prediction. 

The  image  restoration  was  primarily  in  the  3-5  |im  wavelength  region  and  is 
described  in  Appendix  1.  What  is  notable  is  that  the  restoration  work,  based  on 
atmospheric  MTF,  used  prediction  of  turbulence  MTF  instead  of  direct  measurement. 
Furthermore,  it  was  found  that,  contrary  to  common  impressions,  high  spatial 
frequency  filtering  does  not  necessarily  bring  about  improved  image  quality.  It  is 
necessary  to  include  in  the  atmospheic  MTF  modelling  the  unique  shape  of  the  aerosol 
MTF.  It  is  shown  here  that  high  spatial  frequency  filtering  without  considering  the 
unique  MTF  shape  of  such  image  degradation  is  not  nearly  as  effective  as  when  such 
shape  is  considered.  Such  examples  are  included  in  Appendix  1. 

The  key  to  aerosol  MTF  prediction  is  aerosol  size  distribution  prediction.  A 
model  for  relatively  low  humidity  conditions  (our  summer)  has  indeed  been  developed 
now  and  is  included  here  as  Appendix  2.  It  is  relevant  to  most  weather  situations, 
except  for  extreme  ones  such  as  duststorms  or  very  high  humidity  situations.  Efforts 
are  continuing  now  to  extend  it  to  high  humidity  situations  (our  winter)  for  the  next 
report.  The  present  model  indicates  very  high  correlation  of  prediction  with  experiment 
except  for  such  extreme  conditions,  as  shown  in  typical  examples  in  Appendix  2. 


APPENDIX  1 


Restoration  of  thermal  images  distorted  by  the  atmosphere,  using 
predicted  atmospheric  modulation  transfer  function 


D.  Sadot,  G.  Lorman,  R.  Lapardon,  and  N.  S.  Kopeika 
Department  of  Electrical  and  Computer  Engineering 
Ben-Gurion  University  of  the  Negev 
Beer-Sheva,  Israel 


Abstract 

Restoration  of  thermal  images  distorted  by  the  atmosphere,  detected  with  a  focal  plane 
array  (FPA)  Pt-Si  thermal  imaging  system,  is  presented.  The  restoration  method  is  based  upon 
atmospheric  modulation  transfer  function  (MTF)  analysis.  Using  turbulence  and  aerosol  MTF 
prediction  models,  atmospheric  distortions  and  image  degradation  are  modeled.  Restoration 
results  indicate  significant  improvement  in  image  quality.  However,  it  is  critical  to  include  the 
unique  shape  of  aerosol  MTF  when  modeling  atmospheric  MTF  in  order  to  obtain  good 
restoration. 


Introduction 


Thermal  imaging  systems  have  progressed  significantly  in  recent  years.  New  focal  plane 
array  technology  similar  to  ccd  arrays  in  the  visible  range  permits  high  resolution  imaging. 
However,  when  dealing  with  long  distance  imaging  through  the  atmosphere,  system 
performance  Is  often  limited  due  to  the  atmospheric  distortions.  Usually,  the  atmospheric  effect 
in  the  thermal  range  is  referred  to  in  target  acquisition  models^  as  a  constant  attenuation  via 
the  Beer-Lambert  law.  No  treatment  is  considered  for  the  spatial  frequency  dependence  of  the 
atmospheric  effect  i.e.,  the  atmospheric  thermal  imaging  MTF  (primarily  aerosol  MTF).^-^ 
However,  recent  analysis^’^^applying  some  practical  instrumentation  limitations  to  the 
"classical"  theory,^  predicts  significant  atmospheric  MTF  variation  with  spatial  frequency  in  the 
thermal  range  too.  Knowledge  of  atmospheric  MTF  is  essential  to  real-time  high-resolution 
image  restoration.  '  - 

Results  of  measured  atmospheric  MTFs  in  the  long  wave  8-12  jim  thermal  range  with 
low  resolution  imagers  are  shown  elsewhere,^’^  and  are  in  good  agreement  with  theoretical 
predictions.  Here,  turbulence  and  aerosol  MTFs  in  the  3-5  |im  wavelength  atmospheric 
window  with  higher  resolution  imagers  are  estimated  using  both  turbulence^  and  aerosol  MTF 
prediction  models,  with  the  knowledge  of  standard  meteorological  parameters.  Based  on 
atmospheric  MTF,  restoration  of  images  recorded  along  horizontal  atmospheric  paths  is 
presented,  indicating  significant  improvement  iii  image  resolution.  This  demonstrates  the  wide 
variety  of  applications  that  are  achievable  on  the  basis  of  atmospheric  MTF  analysis  in  the 
infrared. 


Experiment 

The  thermal  IR  imaging  system  used  was  a  CCD6003  Platinum  SUicide  Infrared 
Television  Camera.  It  incorporates  a  488  x  512  element  infrared  CCD  (IRCCD)  image  sensor 


designed  to  provide  general-purpose  passive  day/night  imaging  capability.  Sensor  spectral  band 
is  3-5  |im.  The  camera  is  equipped  with  a  100  mm  focal  length  lens  with  a  6.9  (V)  x  9.2  (H) 
degrees  field-of-view.  An  optional  500  mm  focal  length  lens  is  available  too.  Camera 
integration  time  is  1/30  second.  AH  video  signals  were  digitized  via  a  Matrox  frame  grabber 
and  recorded  on  a  personal  computer  for  further  analysis.  Images  were  recorded  through  an 
open  laboratory  window  at  an  altitude  of  about  15  meters  above  the  ground,  horizontally,  for 
image  paths  between  2  to  7  km.  The  receiver  was  located  pn  stable  marble  optical  tables 
located  inside  the  laboratory,  thus  minimizing  significantly  any  possible  mechanical  vibrations. 

During  the  experiment,  measurements  were  taken  of  standard  meteorological  data  - 
ambient  temperature,  relative  humidity,  solar  flux,  wind  speed,  and  particle  size  distribution. 
The  latter  were  measured  over  0.16  to  10  (am  particulate  radius  using  Particle  Measurement 
System  Inc.  (PMS)  instrumentation  .  . 

Simultaneously,  natural  and  urban  ground  scenes  along  the  line-of-sight  were  recorded, 
containing  blur  caused  by  atniospheric  distortions. 

Using  the  predicted  Atmospheric  MTF,  image  restoration  techniques  were  applied  to 
the  distorted  images,  improving  significantly  image  quality. 

Prediction  models 

A  model  for  predicting  turbulence  MTF  already  exists^  {verified  independently  by  U.S. 
Army  Night  Vision  Laboratory).  Also  available  are  IMTURB  or  PROTURB  (U.S.  Army 
Atmospheric  Sciences  Lab.).  With  the  knowledge  of  standard  meteorological  parameters  such 

as  ambient  temperature,  relative  humidity,  wind  speed,  and  solar  flux,  the  turbulence  strength 
parameter  can  be  predicted  quite  accurately.  Using  prediction,  both  short  and  long 

exposure  turbulence  MTFs^  can  be  evaluated. 

In  principle,  the  knowledge  of  aerosol  size  distribution  and  refractive  index,  together 
with  instrumentation  parameters,  can  be  used  to  predict  aerosol  MTF."^  In  Fig.  1,  a  typical 


aerosol  size  distribution  during  the  time  the  experiment  was  carried  out  is  presented.  For 
particulates  of  radii  smaller  than  0.16  |im,  size  distribution  was  extrapolated  according  to  the 
Air  Force  Geophysics  laboratory  model”^,  since  these  could  not  be  measured  with  this  PMS 
instrumentation.  Mie  phase  function  was  calculated  according  to  the  aerosol  size  distribution. 
Fig.  2  presents  aerosol  phase  function  corresponding  to  the  size  distribution  of  Fig.  1.  The 
scattering  and  absorption  coefficients  were  evaluated  too  according  to  Mie  theory.^  The 
medium's  Mie  parameters,  together  with  instrumentation  parameters  including  field-of-view, 
spatial  frequency  bandwidth,  sensor's  dynamic  range,  and  scene's  target  to  background  SNR, 
were  applied  to  the  Practical  Aerosol  MTF  model.*^  Results  of  this  analysis  are  presented  in 
the  following  figures.  First,  the  imaging  system  MTF  is  presented  in  Fig.  3.  In  Tig.  4  the 
"classical"  aerosol  MTF  together  with  the  system  MTF  are  presented  on  a  log-log  scale.  Fig.  5 
presents  the  Specific  Intensity  Function,  which  is  the  Radiative  Transfer  Equation  solution,  but 
limited  according  to  system  MTF.  According  to  Fig.  5  the  ratio  between  the  unscattered  and 
scattered  components  of  light  actually  recorded  in  the  image  is  relatively  small,  unlike  the 
"classical"  model  prediction.^  The  latter  refers  to  the  unscattered  tight  as  an  angular  delta 
function,  thus  predicting  extremely  higher  value  to  the  unscattered-to-scattered  tight  intensity 
ratio.  The  reduction  of  the  ratio  of  the  unscattered  to  scattered  tight  actually  recorded  in  the 
image  is  due  to  the  system's  spatial  frequency  bandwidth  limitation.  This  limitation  causes  a 
severe  blur  of  unscattered  tight,  reducing  this  ratio  to  only  one  order  of  magnitude,  in 
comparison  to  an  angular  delta  function  as  predicted  by  "classical"  theory.  In  Fig.  6  is  an 
example  of  the  Detected  Specific  Intensity  Function  for  particulate  tight  scatter  after’ applying 
aU  the  imaging  system  limitations  (including  finite  field  of  view,  dynamic  range  and  received 
SNR)..  This  represents  an  atmospheric  aerosol  point  spread  function  much  narrower  than  that 
predicted  by  classical  aerosol  MTF  theory  and  is  due  to  the  instrumentation  which  narrows  by 
several  orders  of  magnitude  the  angles  of  tight  scatter  actually  recorded  in  the  image  from  that 
incident  to  the  receiver.^  The  final  result  of  Practical  aerosol  MTF  theory  is  presented  in  Fig.  7 
(continuous  tine).  In  addition,  presented  also  are  the  system  (dashed)  and  turbulence  (dotted) 


MTFs.  Turbulence  MTF  was  predicted  using  the  prediction  modeP  mentioned  above,  with  the 
use  of  meteorological  data  measured  during  the  experiment.  All  three  MTFs  are  multiplied  in 
order  to  obtain  the  overall  system  and  atmospheric  MTF.  This  overall  MTF  for  6.5  km 
distance  is  presented  in  Fig.  8  and  is  the  MTF  used  for  further  image  restoration  analysis  for 
the  objects  located  at  that  range. 


Image  restoration 

Enhancement  techniques  for  thermal  images  are  summarized  elsewhere.  Here,  we  will 
concentrate  on  deblurring  atmospheric  effects. 

Using  the  predicted  thermal  IR  atmospheric  MTF,  several  filtering  techniques  were 
applied  to  some  thermal  images  recorded  by  the  same  imaging  system  in  order  to  restore  the 
blurred  images.  The  most  successive  restorations  were  obtained  by  using  Wiener^i  and 
Improved  Wiener^^  filters.  Other  spatial  frequency  filters,  such  as  inverse  and  pseudo-inverse 
filters,  achieved  less  successful  results.  Fig.  9  presents  an  image  recorded  by  the  IRCCD 
camera  equipped  with  a  100  mm  focal  length  lens.  The  upper  part  of  the  scene  includes 
buildings  at  a  distance  of  1.5  km.  Atmospheric  mrbulence  was  moderate  and  the  atmosphere 
was  relatively  clear  from  particulates  over  the  relatively  short  pathlength.  Thus  it  is  assumed 
that  most  of  the  image  blur  was  caused  by  the  imaging  system  MTF.  In  Fig.  10  is  the  restored 
image  using  the  standard  Wiener  filter  technique.  Though  image  degradation  was  not  severe  as 
explained  above,  there  is  a  distinct  improvement  in  image  quality  after  restoration.  A  good 
indication  of  this  is  the  windows  in  the  remote  buildings  at  the  top  which  are  clearly  sharpened. 

In  general,  if  the  noise  was  negligible,  the  Wiener  filter  would  have  reduced  to  an 
inverse  fdter.  However,  for  any  practical  system,  including  imaging  systems  and  the 
atmosphere,  its  transfer,  function  goes  to  zero  at  high  spatial  frequencies,  so  the  inverse  filter 
will  approach  infinity  at  high  spatial  frequencies.  Since  any  practical  image  contains  noise,  the 
Wiener  filter  limits  image  enhancement  at  high  spatial  frequencies  so  as  to  prevent  infinite 


noise  enhancement  in  the  restored  image.  This  is  demonstrated  in  Fig.  11,  which  represents  the 
standard  Wiener  filters  for  different  SNRs,  such  as  used  in  the  restoration  of  Fig.  9. 

When  atmospheric  conditions  become  worse  or,  conversely,  over  longer  distances, 
atmospheric  limitations  dominate.  Such  an  example  is  demonstrated  in  Fig.  7  where,  for  an 
optical  pathiength  of  6.5  km,  turbulence  and  aerosol  MTFs  are  much  more  severe'than  the 
system  MTF.  An  example  of  a  recorded  image  in  the  IRCCD  camera  equipped  with  a  500 
focal  length  lens,  through  an  optical  pathiength  of  6.5  km,  is  presented  in  Fig.  12.  Predicted 
turbulence  and  aerosol  MTFs  are  very  similar  to  those  presented  in  Fig.  7.  Restoration  of  the 
image  of  Fig.  12  using  the  improved  Wiener  filter  is  presented  in  Fig.  13.  Very  similar  results 
yielding  almost  same  image  quality  were  obtained  using  the  standard  Wiener  filter.  This  is  not 
surprising  since  the  improved  Wiener  filter  is  based  upon  modeling  correctly  the  relative  part 
of  image  degradation  caused  by  turbulence  jitter,  as  explained  in  Ref.  12.  However  here,  due 
to  relatively  low  angular  spatial  frequencies  (up  to  7800  cycles/radian),  and  weak  turbulence 
strength  (about  3-10  [m  ]),  turbulence  jitter  contribution  to  image  degradation  is 
weak.  Therefore,  both  standard  Wiener  and  improved  Wiener  filters  yield  almost  the  same 
results. 

Quantitative  investigation  of  image  correction  is  demonstrated  in  Fig.  14.  The  dashed 
curve  represents  the  overall  (system  and  atmospheric)  predicted  MTF  for  the  image  of  Fig.  12. 
The  dashed-dotted  curve  represents  the  improved  Wiener  filter  used  to  obtain  the  restored 
image  of  Fig.  13.  The  continuous  curve  represents  the  restored  overall  MTF  obtained  after 
image  correction.  Improvement  of  image  quality,  especially  at  high  spatial  frequencies  is 
obvious,  thus  improving  significantly  image  resolution.  This  means  that  essentially  most 
atmospheric  blur  is  removed  in_the  restoration  process,  and  the  distant  object  scene  is  observed 
as  if  there  were  almost  no  atmosphere. 

A  very  important  issue  is  the  restoration  sensitivity  to  blur  estimation  model.  As 
explained  previously,  the  blur  model  includes  three  components;  turbulence  MTF,  practical 
aerosol  MTF,  and  hardware  MTF  (Fig.  7).  Failure  in  modeling  correctly  one  of  the  three 
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components  will  limit  the  quality  of  the  restored  image.  Often  researchers  have  the  impression 
that  most  high  spatial  frequency  filtering  yields  some  image  restoration.  This  is  not  true.  It  is 
critical  to  consider  the  unique  shape  and  cutoff  frequency  of  aerosol  MTF.  Here,  examples  of 
Wiener  filter  restoration  of  the  scene  in  Fig.  9,  using  inaccurate  MTF  estimations,  are 


presented.  MTF  estimations  were  of  a  Gaussian  form,  with  varying  cutoff  spatial  frequencies: 


MTF  =  exp 


(1) 


where  MTF  is  the  Gausian  MTF  estimator,  Q.  is  angular  spatial  frequency  [cycles-rad'i],  and 
Q.C  is  the  estimator  cutoff  frequency.  In  Fig.  15a,  the  overall  blur  MTF  estimation  is  of  a 


Gausian  form  with  a  cutoff  frequency  of  1300  [cycles/radian].  In  Fig.  15b,  a  similar  estimator 
is  used  with  a  cutoff  frequency  of  4300  [cycles/rad].  It  is  obvious  that  in  both  "restorations" 
image  quality  was  reduced  instead  of  being  restored.  Repeating  the  restoration  process  with 
other  Gausian-form  MTF  estimators  of  varying  cutoff  frequencies  yielded  best  restoration 
results  with  a  cutoff  frequency  of  2300  [cycles/rad].  The  restored  image  is  presented  in  Fig.  16. 
Indeed,  the  best  Gausian  estimator  was  most  similar  to  the  actual  predicted  atmospheric  and 
system  MTF  (of  Fig.  8,  using  turbulence  and  aerosol  MTF  prediction  models).  A  comparison 
between  the  Gausian  and  actual  predicted  MTFs  is  presented  in  Fig.  17.  For  the  Gausian 
estimator,  restoration  was  fairly  good  but  not  as  well  as  the  restoration  in  Fig.  13,  which  is 
based  on  the  actual  atmospheric  MTF  prediction  model.  This  comparison  is  also  a  good 
validation  to  the  MTF  prediction  model  since  the  best  MTF  prediction  model  yields  best  image 
restoration.  In  addition,  prediction  and  restoration  procedure  is  implemented  in  real-time, 
saving  a  lot  of  time  for  trial-and-error  estimation  techniques  of  other  image  sharpening 
algorithms  in  order  to  estimate  correctly  the  cutoff  frequency  of  the  blurring  MTF.  Prediction 
can  only  be  possible  by  considering  the  atmospheric  MTF  properly,  including  the  significant 
role  played  by  aerosol  MTF.  Restoration  based  on  MTF  does  not  depend  on  target  shape.  It  is 
a  fundamental  image  correction,  which  can  be  followed  by  other  image  processing  techniques 
as  required.  Knowledge  of  the  atmospheric  MTF  including  aerosol  MTF  permits  faster  image 
correction*  without  the  need  for  many  iterations. 


Conclusions 


Restoration  of  thermal  images  distorted  by  the  atmosphere,  recorded  with  an  IRCCD 
camera,  is  presented.  Restoration  is  based  upon  spatial  frequency  filters  deriving  from  the 
knowledge  of  the  thermal  IR  atmospheric  MTF.  Such  restoration  is  fundamental  and  yields 
image  quality  limited  essentially  by  hardware,  as  if  there  were  almost  no  atmosphere.  This  is  in 
contradiction  to  the  conventional  approach,  which  usually  regards  the  atmospheric  effect  in  the 
IR  range  as  attenuation  only  with  no  spatial  frequency  dependence  considered.  In  principle, 
atmospheric  MTF  can  be  predicted  according  to  weather  forecast,  and  implemented  in  image 
restoration  accordingly.  This  is  important  when  dealing  with  real-time  image  restoration  so 
that  processing  time  is  limited  to  computational  time  only,  with  no  need  for  simultaneous 
atmospheric  measurements.  Knowledge  of  atmospheric  MTF  is  essential  to  good  cost-effective 
system  design  and  is  very  useful  in  image  restoration  without  any  prior  knowledge  of  the  target 


to  be  viewed. 
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Figure  Captions 

Fig.  1.  Atmospheric  particle  size  distribution.  5^  and  Aa  are  scattering  and  absorption 
coefficients  respectively  at  4  |im  wavelength. 

Fig.  2.  Medium's  phase  function. 

Fig.  3.  MTFoftheIRCCD. 

Fig.  4.  Aerosol  ("classical"  model)  together  with  imaging  system  MTF. 

Fig.  5.  Specific  intensity  function  (Radiative  Transfer  solution)  for  the  aerosol  size 
distribution  of  Fig.  1,  based  on  system  MTF  of  Fig.  3. 

Fig.  6.  Detected  (dashed)  and  "classical"  (continuous)  specific -intensity  functions  for  aerosol 
MTF. 

Fig.  7.  Comparison  between  practical  aerosol  (continuous),  turbulence  (dotted),  and  imaging 
system  (dashed)  MTFs. 

Fig.  8.  Overall  atmospheric  and  system  MTF. 

Fig.  9.  Original  picture  recorded  by  the  IRCCD  camera  equipped  with  the  100  mm  * 
focal  length  lens.  The  remote  buildings  are  at  a  distance  of  1.5  km. 

Fig.  10.  Restoration  of  the  original  image  of  Fig.  9. 

Fig.  11.  The  standard  Wiener  Filters  as  used  in  Fig.  10.  From  upper  to  lower  curve  SNR  = 

67,  33,  21,  and  16,  respectively. 

Fig.  12.  Original  picmre  recorded  by  the  IRCCD  camera  equipped  with  the  500  mm  focal 
length  lens.  The  object  is  located  at  a  distance  of  6.5  km. 

Fig.  13.  Restoration  of  the  original  image  of  Fig.  11. 

Fig.  14.  Comparison  between  the  overall  (dashed)  and  restored  (continuous)  MTFs,  and.  the 
■  improved  wiener  filter 

Fig.  15a-b.  Restorations  of  the  original  image  of  Fig.  12,  using  Gausian  blur  model  estimation 
with  different  cutoff  spatial  frequencies;  a)  cutoff  frequency  of  1300  [cycles/md].  b) 

cutoff  frequency  of  4300  [cycles/rad]. 

%  *  ' 


Fig.  16.  Same  as  Fig.  15,  for  a  Gausian  MTF  estimator  with  a  cutoff  frequency  of  2300 
[cycles/rad]. 

Fig.  17.  Comparison  between  atmospheric  MTF  prediction  model  (dashed)  and  best  Gausian 
estimator  (continuous). 


Specific  Intensity  [Watt/m^2/Hz/rad] 


Fig.  5 


Detected  (dashed),  and  classical  (continuous)  Specific  Intensity  functions 


Normalized  MTF  Normalized  MTF 
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Aerosol  (continuous),  system  (dashed),  and  turbulence  (dotted)  MTFs 


Fig.  7 


Overall  MTF  (practical  model) 


Fig.  8 
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Original  picture 


Fig.  9 

Wiener  filter 


Fig.  10 


Fig.  12 
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Fig.  16 


Fig.  17 


APPENDIX  2 


Statistical  model  for  aerosol  size  distribution  parameters  according  to  weather  parameters 
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ABSTRACT 

Predictions  of  atmospheric  transmittance  in  desert  aerosol  environments  using  MODTRAN  code  diverge 
significantly  from  measured  data.  Good  prediction  of  the  desert  particulate  size  distribution  is  required  in 
order  to  predict  atmospheric  scattering  and  absorption  parameters.  It  is  also  essential  to  the  prediction  of 
the  aerosol  atmospheric  modulation  transfer  function  which  is  often  the  dominant  component  of  the  overall 
atmospheric  MTF.^  Recently^  an  effort  to  predict  desert  coarse  aerosol  statistics  but  not  size  distribution 
according  to  simple  weather  parameters  has  been  made.  A  quantitative  analysis  of  the  desert  particulate  size 
distribution  models  was  also  performed^.  In  this  research  the  size  distribution  parameters  measured  by 
optical  counters  are  related  to  weather  parameters.  BCnown  statistical  and  analytical  models  such  as 
MODTRAN  relate  the  size  distribution  parameters  only  to  relative  humidity  for  continental  atmospheres. 
Although  humidity  has  a  significant  role  in  the  prediction  of  aerosol  size  statistics,  other  weather 
parameters  can  also  strongly  influence  the  size  distribution  parameters.  Comparisons  such  as  the  above  can 
be  used  to  predict  under  which  conditions  the  MODTRAN  aerosol  models  have  good  of  poor  accuracy.  It  is 
also  hoped  that  they  will  lead  to  improvements  in  MODTRAN,  improving  the  accuracy  of  the  humidity 
dependence  as  well  as  by  incorporating  other  meteorological  parameters  into  the  MODTRAN  prediction 
models. 


2.  INTRODUCTION 

The  aerosols  that  are  suspended  in  the  atmosphere  are  a  major  factor  that  influence  the  quality  of  images 
that  are  propagating  through  the  atmosphere.  Scattering  and  absorption  of  light  rays  by  the  suspended 
aerosols  contribute  to  the  loss  of  contrast  of  the  received  image.  Forward  light  scattering  in  the  direction  of 
the  image  propagation  also  reduces  the  atmospheric  modulation  transfer  function  (MTF)  and  causes  actual 
image  blur. 

Prediction  of  the  aerosol  size  distribution  can  lead  to  the  prediction  of  the  scattering  and  absorption 
coefficients  and  also  the  scattering  phase  function  of  the  atmospheric  propagating  medium.  These 
properties  characterize  the  extinction  coefficient  of  the  atmosphere  as  well  as  the  aerosol  MTF  of  the 
atmosphere.  Models  such  as  that  suggested  here,  and  extension  of  them,  can  lead  to  prediction  of 
atmospheric  transmission  and  aerosol  derived  blur  according  to  weather. 

3.  EXPERIMENT 

The  particulate  size  distribution  of  the  atmosphere  was  measured  over  the  radius  range  of  0.16  to  10  pm 
using  a  PMS  CSASP-100  optical  particle  counter.  The  optical  counter  was  located  at  a  height  of  25  m 


above  ground.  Integration  time  of  half  an  hour  was  chosen  and  every  half  hour  the  average  particulate  size 
distribution  was  recorded.  Meteorological  parameters  were  also  measured  close  to  the  PMS  location  and 
the  half  hour  average  of  the  temperature,  relative  humidity,  wind  speed,  wind  direction  and  solar  flux  were 
recorded.  The  measurements  were  taken  at  Ben-Gurion  University  of  the  Negev  in  Beer-Sheva,  Israel. 
Data  was  collected  from  June  1991  until  March  1994.  Both  the  PMS  counter  and  the  Campbell  Scientific 
meteorological  station  were  connected  to  a  personal  computer. 

Since  a  goal  of  this  research  is  to  relate  the  aerosol  particulate  size  distribution  parameters  to  simple 
weather  parameters,  the  integration  times  selected  for  counting  and  classifying  particles  according  to  size 
have  to  be  shorter  than  the  meteorological  time  constants.  On  the  other  hand,  a  short  integration  time  can 
lead  to  a  particulate  size  distribution  curve  with  insufficient  data  in  which  case  it  is  hard  to  fit  the  measured 
data  to  a  known  size  distribution  model.  In  our  measurements  an  integration  time  of  half  an  hour  was 
chosen  as  a  compromise  between  the  need  for  a  smooth  size  distribution  curve  and  the  need  to  measure  the 
aerosol  statistics  under  stationary  meteorological  conditions. 

Fig.  1.  describes  the  changes  in  the  standard  deviation  of  the  temperature  along  a  typical  summer  day.  The 
standard  deviation  here  is  calculated  on  the  basis  of  measurement  time  periods  of  half  an  hour.  Fig.  2. 
shows  the  changes  in  the  standard  deviation  of  the  relative  humidity  during  the  same  day.  These  figures 
show  that  the  most  rapid  changes  occur  after  sunrise  and  sunset.  In  our  measurements  the  standard 
deviation  of  the  temperature  did  not  exceeded  2°C  and  that  of  the  relative  humidity  did  not  exceeded  5  %. 
In  most  of  the  cases  the  standard  deviations  were  much  less  than  these  values. 


Fig.  1.  Half  hour  standard  deviation  of  the  temperature  as  a  function  of  time  during  a  typical  summer  day. 
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Fig.  2.  Half  hour  standard  deviation  of  the  relative  humidity  as  a  function  of  time  during  a  typical  summer 
day. 


4.  DATA  ANALYSIS 

The  deviations  in  the  aerosol  size  distribution  are  large.  The  sources  of  the  aerosols  in  the  area  of  our 
measurements  consist  of  a  large  variety  of  aerosol  sources.  The  sources  are  urban  and  rural  aerosol, 
marine,  haze,  industrial  and  desert  dust  aerosols  that  are  transferred  from  two  large  desert  areas  during  dust 
storms 

from  the  Sahara  desert  and  from  the  Saudi  Arabian  desert.*^ 

Fig.  3.  describe  two  extreme  aerosol  size  distribution  that  were  measured  during  the  year  1992.  The 
differences  between  the  two  extreme  cases  is  about  two  orders  of  magnitude.  The  lower  curve  is  not 
continuous  because  particles  at  several  radii  were  not  counted  during  an  integration  period  of  half  an  hour. 
Fig.  4.  shows  the  normalized  standard  deviation  of  the  measured  aerosol  radius  at  each  radius  over  the 
three  year  data  collection  period.  The  normalized  standard  deviation  of  aerosol  concentration  ranges  from 
0.5  at  a  radius  of  0.16  |a,m  to  12  at  a  radius  of  10  fxm.  The  aerosol  size  distribution  consists  of  several 
distributions  that  can  be  modeled  by  superposition  of  lognormal  distributions.  ^  The  standard  deviation  of 
the  particle  count  is  lower  at  small  radii  than  at  large  radii  since  the  low  mode  of  the  aerosol  distribution 
consists  mainly  of  local  aerosols  while  the  large  size  aerosol  population  contains  also  particles  that  are 
transported  from  remote  areas.  Also  the  local  aerosols  are  less  affected  by  changes  in  weather  parameters 
while  the  larger  particles  are  very  sensitive  to  wind  speed  and  to  wind  direction  which,  in  some  cases,  jdelds 
information  about  the  remote  aerosol  source. 

A  sihiple  size  distribution  that  can  characterize  the  aerosol  size  distribution  is  the  Junge  power  law 
distribution.^  The  size  distribution  is  defined  by  only  two  parameters,  A  and  a: 

dnldr  =  Ar-°-  (1)  . 

where  r  is  the  particle  radius  and  n  is  the  particulate  number  concentration. 


The  advantage  of  this  size  distribution  is  that  a  measured  size  distribution  can  be  fitted  to  this  model  by  a 
simple  one  variable  linear  regression. 

In  our  analysis  the  measured  size  distributions  were  fitted  to  the  Junge  power  law  model.  Fig.  5.  describe 
the  average  particulate  size  distribution  and  the  power  law  fit.  All  6000  measurements  of  the  aerosol  size 
distribution  from  May  1992  to  October  1992  were  fitted  to  the  power  law  distribution  and  the  values  of  a 
and  A  were  extracted.  The  average  square  of  the  correlation  coefficient  was  0.95. 

The  power  law  size  distribution  parameters  A  and  a  were  related  to  the  weather  parameters  that  were 
measured  nearby  the  particulate  counter  using  multiple  linear  regression  analysis.  The  summer  model 
obtained  for  A  is  given  by: 

A=aQ+a^rh+aj^rh^  +a^rh^  +a^rk^  (2) 

were  rh  is  relative  humidity  in  percents.  The  coefficients  of  the  model  are  given  in  Table  1. 


Table  1.  Summer  model  for  the  prediction  of  the  parameter  A. 


ao 

ai 

__  _ 

a4 

1.917 

-0.128 

4.61-10-3 

6.85-10-5 

3.73-10-7 

The  model  obtained  for  a  was  the  following: 

o.-h^+b]^rh  +  b2rh^  +b^ws  (3) 

were  is  the  average  wind  speed  in  m-s'^  The  coefficients  for  the  prediction  of  a  are  given  in  Table  2. 


Table  2.  Summer  model  for  the  prediction  of  the  parameter  a. 


bn 

bi 

bi 

2.68 

2.3-10-3 

8.16-10-5 

-2.65-10-2 

The  A  parameter  in  this  model  depends  only  on  the  relative  humidity.  Other  dependences  were  weak  and 
insignificant,  a  depends  also  on  the  average  wind  speed.  Fig.  6  describes  the  measured  a  in  an  ascending 
order  and  the  modeled  a.  From  this  figure  it  is  shown  that  there  is  good  agreement  between  the  measured 
parameter  and  the  modeled  parameter  for  about  70%  of  the  measurements.  Fig.  7.  shows  the  measured  and 
modeled  parameters  A  in  ascending  order. 

Figs.  8.  show  some  random  measured  particulate  aerosol  distributions  (solid  lines)  and  the  predicted  size 
distribution  according  to  the  measured  weather  parameters  (dashed  lines).  The  curves  shows  that  there  is 
very  good  agreement  between  the  measured  data  and  the  model  for  the  parameters  a  and  A  for  various 
relative  humidity  conditions.  The  straight  line  represents  results  (1)  -  (3)  and  are  not  a  "best  fit"  to  the 
measured  aerosol  size  distribution. 


Fig.  3.  Average  aerosol  size  distribution  curve  and  size  distributions  for  minimum  and  maximum  numbers 
of  particles  that  were  measured. 


Fig.  4.  Standard  deviation  of  all  aerosol  concentrations  normalized  according  to  mean  aerosol 
concentration.. 


Fig.  5.  Measured  average  particulate  size  distribution  and  power  law  model. 


Fig.  6.  Comparison  of  measured  and  predicted  values  of  a  in  ascending  order. 


Fig.  7.  Comparison  of  measured  and  predicted  values  of  A  in  ascending  order. 


Fig.  8a.  Measured  and  modeled  (straight  line)  size  distributions.  Wd  and  Sf  are  wind  direction  and  solar 
flux,  respectively. 


statistics  of  the  aerosols  are  influenced  strongly  by  the  local  weather.  There  are  few  exceptions,  mainly 
when  the  aerosol  size  distribution  consist  of  transported  particles  such  as  desert  dust  particles  or  marine 
particles  that  are  brought  from  other  areas. 

For  further  investigation,  the  cases  when  the  size  distribution  curve  is  affected  by  remote  particles  should 
be  studied  and  the  sources  of  these  aerosols  should  be  identified  by  analysis  of  the  history  of  the  wind  speed 
and  wind  direction. 
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